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Abstract 

We prove the existence of a large class of asymptotically flat initial data 
with non- vanishing mass and angular momentum for which the metric and 
the extrinsic curvature have asymptotic expansions at space-like infinity 
in terms of powers of a radial coordinate. 
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1 Introduction 



An initial data set for the Einstein vacuum equations is given by a triple 
(S, h a b, ^ab) where S is a connected 3-dimensional manifold, h a b a (positive 
definite) Riemannian metric, and a b & symmetric tensor field on S. The data 
will be called "asymptotically flat" , if the complement of a compact set in S can 
be mapped by a coordinate system x 3 diffeomorphically onto the complement 
of a closed ball in R 3 such that we have in these coordinates 

^• = (l + ^)5« + 0(f- 2 ), (1) 



= 0(f- 2 ), (2) 

as f = (X)j=i(^"') 2 ) 1 ^ 2 ~~ * 00 • Here the constant m denotes the mass of the 
data, a,b,c... denote abstract indices, i,j,k..., which take values 1,2,3, denote 
coordinates indices while Sij denotes the flat metric with respect to the given 
coordinate system x 3 . Tensor indices will be moved with the metric h a b and its 
inverse h ab . We set Xi = x l and d 1 = di. Our conditions guarantee that the 
mass, the momentum, and the angular momentum of the initial data set are 
well defined. 

There exist weaker notions of asymptotic flatness (cf. but they are 

not useful for our present purpose. In this article we show the existence of a 
class of asymptotically flat initial data which have a more controlled asymptotic 
behavior than ([l]), (||) in the sense that they admit near space- like infinity 
asymptotic expansions of the form 

(3) 

k>2 



9i 



k>2 



(4) 



where and ^ are smooth function on the unit 2-sphere (thought as being 
pulled back to the spheres f = const, under the map x? — > x 3 /r). 

We are interested in such data for two reasons. In pj| has been studied in 
considerable detail the evolution of asymptotically flat initial data near space- 
like and null infinity. There has been derived in particular a certain "regularity 
condition" on the data near space-like infinity, which is expected to provide a 
criterion for the existence of a smooth asymptotic structure at null infinity. To 
simplify the lengthy calculations, the data considered in |^3| have been assumed 
to be time-symmetric and to admit a smooth conformal compactification. With 
these assumptions the regularity condition is given by a surprisingly succinct 
expression. With the present work we want to provide data which will allow 
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us to perform the analysis of |23|] without the assumption of time symmetry 
but which are still "simple" enough to simplify the work of generalizing the 
regularity condition to the case of non-trivial second fundamental form. 

Thus we will insist in the present paper on the smooth conformal compact- 
ification of the metric but drop the time symmetry. A subsequent article will 
be devoted to the analysis of a class of more general data which will include in 
particular stationary asymptotically flat data. 

The "regular finite initial value problem near space-like infinity" , formulated 
and analyzed in Q , suggests to calculate numerically entire asymptotically flat 
solutions to Einstein's vacuum field equations on finite grids. In the present 
article we provide data for such numerical calculations which should allows us 
to study interesting situations while keeping a certain simplicity in the handling 
of the initial data. 

The difficulty of constructing data with the asymptotic behavior (||), (||) 
arises from the fact that the fields need to satisfy the constraint equations 

D b ^ ab -D a V = 



R + ^ 2 - * afc * ab = 0, 

on S, where D a is the covariant derivative, R is the trace of the corresponding 
Ricci tensor, and ^ = h ah ^ a b- Part of the data, the "free data", can be given 
such that ([|), (^) hold. However, the remaining data are governed by elliptic 
equations and we have to show that (||), ([!]) are in fact a consequence of the 
equations and the way the free data have been prescribed. 

To employ the known techniques to provide solutions to the constraints, we 
assume 

* = 0, (5) 

such that the data correspond to a hypersurface which is maximal in the solution 
space-time. 

We give an outline of our results. Because of the applications indicated 
above, we wish to control in detail the conformal structure of the data near 
space-like infinity. Therefore we shall analyze the data in terms of the conformal 
compactification (S,h a b,^ a b) of the "physical" asymptotically flat data. Here 
S denotes a smooth, connected, orientable, compact 3-manifold. It contains a 
point i such that we can write S = S\{i}. The point i will represent, in a sense 
described in detail below, space like infinity for the physical initial data. 

By singling out more points in S and by treating the fields near these points 
in the same way as near i we could construct data with several asymptotically 
flat ends, since all the following arguments equally apply to such situations. 
However, for convenience we restrict ourselves to the case of a single asymptot- 
ically flat end. 

We assume that h a b is a positive definite metric on S with covariant deriva- 
tive D a and ^ a b is a symmetric tensor field which is smooth on S. In agreement 
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with (|^) we shall assume that ^ ab is trace free, 

h ab V ab = 0. 

The fields above are related to the physical fields by rescaling 

hab = # 4 h ab , * ab = 9- 2 * ab , (6) 

with a conformal factor 9 which is positive on S. For the physical fields to 
satisfy the vacuum constraints we need to assume that 

D a $ ab = on S, (7) 

(D b D b -lR)9 = -i* ab ^ b 9- 7 on S. (8) 
o o 

Equation (||) for the conformal factor 9 is the Lichnerowicz equation, transferred 
to our context. 

Let x J be /i-normal coordinates centered at i such that hki = Ski at i and set 
r = {J^i^ii^) 2 ) 1 ^ 2 ■ To ensure asymptotical flatness of the data ^ we require 

^ ab = 0(r- 4 ) as r^O, (9) 
limr0 = l. (10) 

In the coordinates x j = /r 2 the fields (§) will then satisfy (0), (|) (cf. ||, 
[ p3| for this procedure). 

Not all data (0) derived from data h ab , ^ a ft as described above will satisfy 
conditions (^), (^. We will have to impose extra conditions and we want to 
keep these conditions as simple as possible. 

Since we assume the metric h ab to be smooth on S, it will only depend on 
the behavior of 9 near i whether condition (^|) will be satisfied. Via equation (||) 
this behavior depends on ^ab- What kind of condition do we have to impose on 
"fab in order to achieve (||) ? 

The following space of functions will play an important role in our discussion. 
Denote by B a the open ball with center i and radius r — a > 0, where a is chosen 
small enough such that B a is a convex normal neighborhood of i. A function 
/ G C°°(>S) is said to be in E°°(B a ) if on B a we can write / = A + rf 2 with 



/i,/ 2 € C°°{B a ) (cf. definition 3.7) 



An answer to our question is given by following theorem. 

Theorem 1.1 Let h ab be a smooth metric on S with positive Ricci scalar R. 
Assume that ty ab is smooth in S and satisfies on B a 

r 8 * ab * ab EE°°{B a ). (11) 

Then there exists on S a unique solution 9 of equation (j^), which is positive, 
satisfies and has in B a the form 

6 = -, 9eE°°(B a ), 9(i) = l. (12) 
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In fact, we will get slightly more detailed information. We find that 6 = 
Ui +ru2 on B a with U2 E E°°(B a ) and a function u\ E C°°(B a ) which satisfies 
til = 1 + 0{r 2 ) and 



-.6 1 „s u l 



(D b D — -R) — 



r 



1R, 



in B a \{i}, where 9r is in C co (Ba) and vanishes at any order at i. 

If 9 has the form (jl^) then (||) will be satisfied due to our assumptions on 

h a h. 

Note the simplicity of condition (|IT|). To allow for later generalizations, we 
shall discuss below the existence of the solution 9 under weaker assumptions on 
the smoothness of the metric h a b and the smoothness and asymptotic behavior 



of (cf. theorem 3.4). In fact, already the methods used in this article would 
allow us to deduce analogues of all our results under weaker differentiability 
assumptions, however, we are particularly interested in the C°° case because it 
will be convenient in our intended applications. If the metric is analytic on B a 
it can be arranged that 9r = and u\ is analytic on B a (and unique with this 
property, see |^4| and the remark after Theorem |3.12|) . We finally note that the 
requirement R > 0, which ensures the solvability of the Lichnerowicz equation, 
could be reformulated in terms of a condition on the Yamabe number (cf . J2?J ) . 
It remains to be shown that condition ( |Tl| ) can be satisfied by tensor fields 
which satisfy (|^), @. A special class of such solutions, namely those which 
extend smoothly to all of S, can easily be obtained by known techniques (cf. 
[fl6|). However, in that case the initial data will have vanishing momentum 
and angular momentum. To obtain data without this restriction, we have to 
consider fields ^ a b E C°°(S) which are singular at i in the sense that they admit, 
in accordance with (||), (^), (jTo|), at i = {r = 0} asymptotic expansions of the 
form 

*y ~ J2 ^ k if k with *y G C°°{S 2 ). (13) 

k>-4 

It turns out that condition (|ll| ) excludes data with non- vanishing linear mo- 
mentum, which requires a non- vanishing leading order term of the form 0(r~ 4 ). 



In section 3.4 we will show that such terms imply terms of the form log r in 



9 and thus do not admit expansion of the form (H). However, this does not 



necessarily indicate that condition (11) is overly restrictive. In the case where 



the metric h a b is smooth it will be shown in section 3.4 that a non-vanishing 
linear momentum always comes with logarithmic terms, irrespective of whether 
condition JTl| ) is imposed or not. 

There remains the question whether there exist fields ^ a b which satisfy ( p"T| ) 
and have non-trivial angular momentum. The latter requires a term of the form 



0(r 3 ) in ([L3ft. It turns out that condition (11) fixes this term to be of the form 



A 3 
^f/ = ■^(3n i n j ~ <%) + - s (n j t ka J l n k + ni ei jk J k n l ), (14) 



5 



where n l = x l /r is the radial unit normal vector field near i and J k , A are 
constants, the three constants J k specifying the angular momentum of the data. 
The spherically symmetric tensor which appears here with the factor A agrees 
with the extrinsic curvature for a maximal (non-time symmetric) slice in the 
Schwarzschild solution (see for example (T(]] ) . Note that the tensor ^fj J satisfies 

condition ( pd| ) and the equation d 1 ^^' 1 — on S for the flat metric. In the next 
theorem we prove an analogous result for general smooth metrics. 

Theorem 1.2 Let h a b a smooth metric in S. There exists trace- free tensor 
fields 'J'ab G C°°(S\ {i}) satisfying §Ufy with the following properties 

(i) * afc = *f b J + * ab , where tf^ls given by @) and * ab = 0(r- 2 ). 

(ii) D a m ab = onS 

(Hi) r 8 ' 1 $> a i ) 1 $> ab satisfies condition ftH\). 



We prove a more detailed version of this theorem in section (4.3). There it 



will be shown how to construct such solutions from free-data by using the York 



splitting technique (|p5|). In section 4T the case where h a b is conformal to the 



Euclidean metric is studied in all generality. 



2 Preliminaries 

In this section we collect some known facts from functional analysis and the 
theory of linear elliptic partial differential equations. 

Let Z be the set of integer numbers and No the set of non negative integers. 
We use multi-indices (3 = (ft, ft, . ..,ft) G N% and set |ft = £™ =1 ft, ft = 
ft!ft!...ft!, = (x 1 )' 3l (x 2 ) fe ...(x n )^, d^u = d^d^ ...d^u, D^u = 
... Dfr u, and, for ft 7 G N», /? + 7 = (ft + 7l , . . . , ft + 7n ) and (3 < 7 
if ft < 7i< We denote by SI an open domain in K 3 (resp. in S; quite often we 
will then choose SI = B a ). 

We shall use the following functions spaces (see jj), (^H) for definitions, no- 
tations, and results): the set of m times continously differentiable functions 
C m (S7), the Holder space C m,a (S7), where < a < 1, the corresponding spaces 
C" l (Sl), C m,a (Sl), the space Cq 30 (SI) of smooth function with compact support in 
S7, the Lebesgue space L p (£l), the Sobolev space W m ' p (£l), and the local Sobolev 
space W^' p (f2). For a compact manifold S we can also define analogous spaces 
LP(S), C m < a {S), W m 'P(S) (cf. §). 

We shall need the following relations between these spaces. 

Theorem 2.1 (Sobolev imbedding) Let £1 be a C ' 1 domain in K 3 , let k, 
m, j be non-negative integers and 1 < p, q < oo. Then there exist the following 
imbeddings: 

(i) Lf mp < 3, then 

w j+m, P ^ c W j, q ^ p<q< 3p/ (3 - mp) . 
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(ii) If (m — l)p < 3 < mp, then 



W 3+ m , P (n) c C j> (o) ) a = m - 3/ p . 



Theorem 2.2 Let u G W ' (fl), and suppose there exist positive constants a < 
1 and K such that 

I \du\dfi< KR 2+a for all balls B R C ft o/ radius R > 0. 
TTien u e C Q (ft). 

Our existences proof for the non-linear equations relies on the following 
version of the compact imbedding for compact manifolds [^| . 

Theorem 2.3 (Rellich-Kondrakov) The following imbeddings are compact 

(i) W m ' q {S) C D>(S) if I > l/p> l/q-m/3 > 

(ii) W m - q {S) c C a {S) ifm-a> 3/q, < a < 1. 

A further essential tool for the existence proof is the Schauder fixed point the- 
orem [ p5| . 

Theorem 2.4 (Schauder fixed point) Let B be a closed convex set in a Ba- 
nach space V and let T be a continuous mapping of B into itself such that the 
image T(B) is precompact, i.e. has compact closure in B. Then T has a fixed 
point. 

We turn now to the theory of elliptic partial differential equations (see p"2| , 
[ fl4| , p5[ , pl|). Let L be a linear differential operator of order m on the compact 
manifold S which acts on tensors fields u. In the case where u ~ u° 1 " ,a " , » is a 
contravariant tensor field of rank mi, L has in local coordinates the form 

m m 

Lu= ]T a^-^ ll ... lmi0 D^u^-^ = ]T apD^u, (15) 

1(91=0 |0|=O 

where the coefficients aP 1 '"^' i,... !mi ^ = a(s) J:L " Jm 2 ^..j p are tensor fields of 
a certain smoothness, and D denotes the Levi-Civita connexion with respect to 
the metric h. In the expression on the right-hand side we suppressed the indices 
belonging to the unknown and the target space. Assuming the same coordinates 
as above, we write for given covector at a point x G ft and multi-index (3 as 
usual £0 = £f 1 . . . £f 4 and define a linear map A(x, £) : R™ 1 -> R™ 2 by setting 
(i(i,0«) J1 '' J '" 2 = E^i^aW 1 -^ 2 4l ... Imi /3^u ll - 4mi - The operator L is 
elliptic at x if for any £ ^ the map A(x,£) is an isomorphism, L is elliptic on 
S if it is elliptic at all points of S. 

We have the following W regularity result §, §, Q, |l). 
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Theorem 2.5 (L p regularity) Let L be an elliptic operator of order m on fl 
(resp. S) with coefficients ap £ W S ^^' P (S), where > 3/p + k — m + 1, and 
p > 1. Let s be a natural number such that Sk > s — m > 0. Let u £ 1^;™ C ' P (^) 
(resp. WT?' P (S) ), with p > 1, be a solution of the elliptic equation Lu = /. 

(i) Iffe Wj*- m ' 9 (n), q > p, then u £ W£«(Q). 

(nj I// E W s ' m ' q (S), q>p, then u £ W s ' q (S). 

Furthermore, we have the Schauder interior elliptic regularity [|| , JlS| ] , p5[ , 

[ED- 

Theorem 2.6 (Schauder elliptic regularity) Let L 6e an elliptic operator 
of order m on U with coefficients ap £ C k,a (pl). Let u £ W m,v (VL), with p > 1, 
be a solution of the elliptic equation Lm = /, with f £ C > a (£l). Then u £ 
C k+m < a {VL'), for all fi' CC Q. 

For linear elliptic equations we have the Fredholm alternative for elliptic 
operators on compact manifolds Il2j . 

Theorem 2.7 (Fredholm alternative) Let L be an elliptic operator of order 
m on S whose coefficients satisfy the hypothesis of Theorem \2.5\ Let s be some 
natural number such Sfc > s — m > and f £ L P (S), p > 1. Then the equation 
Lu = f has a solution u £ W m ' p (S) iff 

/ < v, f >h dfj, = for all v £ ker(L*). 
Js 

Here dfi denotes the volume element determined by h and L* the formal adjoint 
of L, which for the operator ( |l5| ) is given by 

m 

L*u= (-l) m D^(a u). (16) 

|/3|=0 

Furthermore <,>h denotes the appropriate inner product induced by the metric 
h a b. In our case, where u and / will be vector fields /° and u a , we have 
<uj > h = f a u a . 
Let 

Lu = diicPdjU + b l u) + dd l u + du, (17) 

be a linear elliptic operator of second order with principal part in divergence 
form on f2 which acts on scalar functions. An operator of the form (|l7j ) may be 
written in the form ( |l5|) provided its principal coefficients a y are differentiable. 

We shall assume that L is strictly elliptic in il; that is, there exists A > 
such that 

a ij (x)^j > A|£| 2 , VseO, eeR™. (18) 
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We also assume that L has bounded coefficients; that is for some constants 
A and v>0we have for all x E ft 

^|a y '| 2 <A 2 , A- 2 5^(|6 i | 2 + |c i | 2 ) + A- 1 |d|<iy a . (19) 

In order to formulate the maximum principle, we have to impose that the 
coefficient of u satisfy the non-positivity condition 

(dv - WdiV) dx < Vt)>0,D6Cj(O). (20) 



n 

We have the following versions of the maximum principle p5[ . 

Theorem 2.8 (Weak Maximum Principle) Assume that L given by \l r , 
satisfies conditions filfy, Q ) and Let u e V^ 1,2 (S1) satisfy Lu > (< 0) 

in tt. Then 

sup u < sup u + (inf u > inf it - ) , 
n on n an 

where u + (x) = max{u(a;), 0}, vr{x) — min{w(x),0}. 

Theorem 2.9 (Strong Maximum Principle) Assume that L given by J? 
satisfies conditions ([7^j, and (p^j. Let it S W 1,2 (il) satisfy Lu > in fl 

Then, if for some ball B CC f2 we /iaz;e 

sup w = supu > 0, 

B fi 
i/ie function u must be constant in fl. 

Because u is assumed to be only in W 1 ' 2 the inequality Lu > has to be 
understood in the weak sense (see [f25| for details). 

3 The Hamiltonian Constraint 



In this section we will proof theorem (1.1) 



3.1 Existence 

The existence of solutions to the Lichnerowicz equation has been studied under 
various assumptions (cf. [|l5), juSI, (28) and the reference given there). The 
setting outlined above, where we have to solve (||), ( |l0| ) on the compact manifold 
S, has been studied in §, @, Q. 

In general the "physical" metric provided by an asymptotically flat initial 
data set will not admit a smooth conformal compactification at space-like infin- 
ity. Explicit examples for such situations can be obtained by studying space-like 
slices of stationary solutions like the Kerr solution. To allow for later general- 
izations of the present work which would admit also stationary solutions we 
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shall prove the existence result of theorem 1.1 for metrics h a b which are not 
necessarily smooth. The proof we will employ Sobolev spaces W m ' p (S) and 
the corresponding imbeddings and elliptic estimates (in particular, there will be 
no need for us to employ weighted Sobolev spaces with weights involving the 
distance to the point i). With these spaces and standard LP elliptic theory we 
will also be able to handle the mild r _1 -type singularity at i which occurs on 
the right hand side of equation (||). 

The conformal Laplacian or Yamabe operator 

L h = h ab D a D b - Ir, 
o 

which appears of the left hand side of d|), is a linear elliptic operator of second 
order whose coefficients depend on the derivatives of the metric h up to second 
order. The smoothness to be required of the metric h is determined by the 
following considerations. In the existence proof we need: 

(i) The existence of normal coordinates. This suggests to assume h S C 1 ' 1 (S I ). 



(ii) The maximum principle, theorems ^8 and |2.9| . The required boundedness 



of the Ricci scalar R imposes restrictions on the second derivative of h. 



(iii) The elliptic LP estimate, theorem 2J5. This requires that h £ W 3 ' P (S) for 
p > 3/2. 

Since the right-hand side of equation (@) is in L? the assumption that h 6 
W 3 ' P (S), p > 3, would be sufficient to handle equation (||). However, when we 



will discuss the momentum constraint in section 4.2 , we will wish to be able 
to handle cases where p < 3/2. In these cases the conditions of theorem |2.5| 
suggest to assume that h £ W 4,P (S). In order to simplify our hypothesis we 
shall assume in the following that 

h ab e W 4 *(S), p>3/2. (21) 

The imbedding theorems then imply that h a b G C 2 ' a (S), < a < 1, whence 
R G C a (S). We note that ( ^lf ) is not the weakest possible assumption but it 
will be sufficient for our future applications. 



Lemma 3.1 Assume that h satisfies l[2\ ) and R > 0. Then: 

(i) Lh : W 2,q (S) — > L q (S), q > 1, defines an isomorphism. 

(ii) if u G W 1,2 (S) and L^u < 0, then u > 0; if, moreover, L^u / £ L q (S), 
then u > 0. 

Proof: (i) To show injectivity, assume that LhU = 0. By elliptic regularity u 
is smooth enough such that we can multiply this equation with u and integrate 
by parts to obtain 



j (r> a uD a u + Iru 2 ^) dHh - o. 
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Since R > it follows that u — 0. Surjectivity follows then by theorem |2.7| since 
Lh = L* h . Boundedness of Lh is immediately implied by the assumptions while 
the inequality 

IMIw 2 .»(S) < C\\L h u\\ LP{s) , 



which follows from the elliptic estimates underlying theorem 2.5 and the injec 



tivity of Lh (see e.g JL21 for this well known result), implies the boundedness of 
L h 1 - 

(ii) If we have u < in some region of S, it follows that sup s (— u) > 0. 
Then there is a region in S in which we can apply the maximum principle to 
the function —u to conclude that u must be a non-positive constant whence 
LhU = —Ru/8 > in that region. In the case where L^u < we would arrive 
at a contradiction. In the case where L^u < we conclude that u — in the 
given region and a repetition of the argument gives the desired result. H 



To construct an approximate solution we choose normal coordinates x 3 cen- 
tered at i such that (after a suitable choice of a > 0) we have in the open ball 
B a in these coordinates 

hij = S i:j + hij, h ij = 5 ij + h ij (22) 

with 

hij = 0(r 2 ), h ij = 0(r 2 ), x%j = 0, Xl h ij = 0. 

Notice that h%j, h K \ defined by the equations above, are not necessarily related 
to each other by the usual process of raising indices. 

Denoting by A the flat Laplacian with respect to the coordinates x 3 , we 
write on B a 

L h = A + L h , 

with 

L h =h ij d i dj+b i d i -\R. (23) 
o 

We note that 

b l =0{r). 

Choose a function \a € C°°(S) which is non- negative and such that \a = 1 
in B a / 2 and Xa = in S \ B a . Denote by <Jj the Dirac delta distribution with 
source at i. 

Lemma 3.2 Assume that h satisfies and R > 0. Then, there exists a 
unique solution 0q of the equation L^Oq = — 47r<5,;. Moreover 9q > in S and we 
can write 0q — Xa/r + g with g 6 C a (S), < a < 1. 
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Proof: Observing that 1/r defines a fundamental solution to the flat Laplacian, 
we obtain 

A(^) = -47r5 i+ x, (24) 
r 

where x is a smooth function on S with support in B a \ B a / 2 - The ansatz 
$o = Xa/f + 9 translates the original equation into an equation for g 

L h g = -L h {—) - x- 
r 



A direct calculation shows that Lh{Xa r ^ r ) £ L q (S), q < 3. By lemma B.l there 
exists a unique solution g £ W 2,q (S) to this equation which by the imbedding 
theorem is in C a (S). 

To show that 8q is strictly positive, we observe that it is positive near i (be- 
cause r _1 is positive and g is bounded) and apply the strong maximum principle 
to -0 O . ■ 

We use the conformal covariance of the equation to strengthen the result on 
the differentiability of the function g. Consider a conformal factor 

luo = e /o with f £ C°°(S) such that / = - x 3 x k L jk (i) on B a , (25) 

where we use the normal coordinates x k and the value of the tensor 

Lab = Rab — -^Rh a b, (26) 

at i. Then the Ricci tensor of the metric 

Kb = 4hab (27) 

vanishes at the point i and, since we are in three dimensions, the Riemann 
tensor vanishes there too. Hence the connection and metric coefficients satisfy 
in the coordinates x k 

T^" fc = 0(r 2 ), h' v] , = Sij + 0(r 3 ). (28) 



Corollary 3.3 The function g found in lemma 3.i is in C 1 ' a (S), < a < 1. 
Proof: With uq = e? and h' ab as above, we note that 

L h ,(6' )=Uo 5 L h (6 ), (29) 

where 

O'o = u>o%. (30) 
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We apply now the argument of the proof of lemma |3.2| to the function 9' . 
Since we have by equation (p8|) that Lh'ixar^ 1 ) € L°°(S), it follows that 



'o 



Xa/r + g', where g' e C ,a (S). We use equation (30), the fact that 



loq = 1 + 0(r 2 ), and lemma |3^ to obtain the desired result. 
We note that the function 



0-1 
7 o 



Xa + rg 



(31) 



is in C a (S), it is non-negative and vanishes only at i. To obtain 9, we write 
9 = 9q + u and solve on S the following equation for u 

L h u = -~^ 7 * ai) * ob (l + 9^u)- 7 . (32) 



Theorem 3.4 Assume that h ab e W 4 ' P (S) with p > 3/2, that R > on S, 
and that 9Q 7 ^ ab ^! ah S L q (S), q > 2. T/ien i/iere exists a unique non-negative 
solution u £ W 2 ' q (S) of equation We We ti > on S un/ess * a 6 v I /a '' = 

6l«(S'). 

We note that our assumptions on ^> a i, impose rather mild restrictions, which 
are, in particular, compatible with the fall off requirement (m. By the imbedding 



theorem gj] we will have u € C a (S), a = 2 - 3/q, for 9 > 3; and u € C 1 '"^), 
for g > 3. 

Proof: The proof is similar to that given in ||, with the difference that we 
impose weaker smoothness requirements. Making use of lemma [O], we define 
a non-linear operator T : B — > C° (S 1 ) , with a subset _B of C° (S 1 ) which will be 
specified below, by setting 

T(u)=L^f(x,u), 

where 

f(x 1 u) = -^ 9 7 ^ ab ^ ab g(x 1 u) (33) 

with 

sex, u) = (l + ^r 7 . 

In the following we will suppress the dependence of / and g on x. Let if) € 
W 2 ^(S) C C a {S) be the function satisfying ijj = T(0) and set B = {u € 
C°(S) : < u < ip}, which is clearly a closed, convex subset of the Banach 
space C°(S). 

We want to use the Schauder theorem to show the existence of a point u E B 
satisfying u — T(u). This will be the solution to our equation. 
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We show that T is continuous. Observing the properties of B^ 1 noted above, 
we see that g defines a continuous map g : B — > L 2 . Using the Cauchy-Schwarz 
inequality, we get 

ll/K) - f(u 2 )\\ L , < ||i^ 7 * ab * ab || L2 || 5 K) - g{u2)\W- 



By theorems |j, |j we know that the map L 2 (S) -»■ W 2 ' 2 (S) -> C ^), where 
the first arrow denotes the map L^ 1 and the second arrow the natural injection, 
is continuous. Together these observations give the desired result. 

We show that T maps B into itself. If u > we have / < whence, by 



lemma 3.1, T{u) = L~\f{u)) > 0. If u x > u 2 it follows that /(tt x ) > f{u 2 ) 



whence, again by lemma 3.1, T(u 2 ) — T{u\) = L h {f{u 2 ) — /(tti)) > 0. We 



conclude from this that for u e B we have < < T(0) = ip. 

Finally, T(B) is p recompact because W 2 ' 2 (S) is co mpa ctly embedded in 



C°(S) by theorem ^3. Thus the hypotheses of theorem ( |2.4|) are satisfied and 
there exists a fix point it of T in _B. By its construction we have u G L q (S) 
whence, by elliptic regularity, u £ W 2 ' q (S). 

To show its uniqueness, assume that Ui and u 2 are solutions to (|3^). Ob- 
serving the special structure of g and the identity 

6 

a' 7 -c- 7 = {c-a) JV-V 1 -*, 
i=o 

which holds for positive numbers a and c, we find that we can write Lh(ux— u 2 ) — 
c(ui — u 2 ) with some function c > 0. The maximum principle thus allows us 
to conclude that u\ — u 2 . The last statement about the positivity of u follows 



from lemma 3.1, (ii). 



3.2 Asymptotic expansions near i of solutions to the Lich- 
nerowicz equation 

The aim of this section is to introduce of the functional spaces E m ' a , to point 



out the simple consequences listed in Lemma |3.9|, and to proof Theorem 3.11 



These are the tools needed to prove Theorem 1.1 



Let m £ No, and denote by V m the space of homogenous polynomials of 
degree m in the variables xK The element of V m are of the form X)|/3|=m x ^ 
with constant coefficients Cp. Note that r 2 is in V 2 but r is not in V\. We 
denote by 7i m the set of homogeneous harmonic polynomials of degree to, i.e 
the set of p € V m such that Ap = 0. For s € Z, we define the vector space r s V m 
as the set of functions of the form r s p with p € V m . 

Lemma 3.5 Assume s€Z. The Laplacian defines a bijective linear map 

A : r s V m -» r s - 2 V m , 

in either of the following cases 

(i) s > 

(ii) s < 0, \s\ is odd and to + s > 0. 
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Note that the assumptions on m and s imply that the function A(r s p m ) G 
C°°(R 3 \ {0}) defines a function in Lj oc (R 3 ) which represents A(r s p m ) in the 
distributional sense. 

Proof: Since A maps V m into V m -2, we find from 

A(r s p) = r s - 2 (s(s + 1 + 2m)p + r 2 Ap) , (34) 

that A(r s p) G r s - 2 V m - 

We show now that the map is bijective for certain values of s and m. Because 
r s V m and r s ~ 2 V m have the same finite dimension, we need only show that the 
kernel is trivial for some s and m. The vector space V m can be written as a 
direct sum 

V m =H m (5 r 2 H m - 2 © r 4 H m -4 ■ • • , (35) 
(cf. |(J). If A(r s p) = 0, we get from (|§) that 

= A(r s p) = J2 &(r s+2k h m -2k), 

0<k<m/2 

with h m -2k € H m -2k- Applying (g), we obtain 

0= ^ r 2fc (s + 2fc)(s + l + 2(m-fc))/ lm „ 2fc , 

0<fc<m/2 

which allows us to conclude by ([35]) that 

(s + 2k)(s + 1 + 2(m - k))h m -2k = 0. 

Since by our assumptions (s + 2fc)(s + 1 + 2(m — fc)) 7^ 0, it follows that the 
polynomials h m —2k vanish, whence r s p = 0. I 

We will need the following technical Lemma regarding Holder functions: 

Lemma 3.6 Suppose m G N, < a < 1, / G C m ' a (B a ), and T m denotes the 
Taylor polynomial of order m of f. Then f R = f — T m is in C m ' a (B a ) and 
satisfies, if \(3\ < m, 

d f R = 0(r m -^ l+a ) as r^0. 

Moreover, let s be an integer such that s < 1 and m + s — 1 > 0. TTien /r 
satisfies: 

(i) r s ~ 2 f R G W m+a - x *(B e ), /orp< 3/(1 - a), < e < a. 

(ii) Ifm + s-l>l then r s ~ 2 f R G C m+s_2,a (B e ). 
fit»; r/n G C m ' a (B e ). 

Proof: The relation 

\d p f R \ < C\x\ m -W +a , x G S e , (36) 



15 



is a consequence of Lemma 

(i) We have 

d^r s - 2 f R )= 0^'f R d^'(r s - 2 ), 
with certain constants Cpi , and the derivatives of r s ~ 2 are bounded for ieB e 

by 

|a 7 'r s - 2 | < Cr s - 2 - |Vl . 

Observing (|36|), we obtain 

for ie5 e , whence, by our hypothesis m + s — 1 > 0, 

\d P {r s - 2 f R )\ <Cr- 1+a , (37) 

for < m+s—1. Using that r~ 1+Q is in L p (B e ) forp < 3/(1 — a), we conclude 
that d (r s ~ 2 f) £ L p (B e ) whence r^ 2 /^, € W m+s ~ 1,p (B e ). 

(ii) From the relation above and theorem 2.1 we conclude for m + s — 1 > 1 
that r s ~ 2 / fl e C m+s - 2 ' a< '(B e ) with a' = 1 - 3/p < a. To show that a 1 can 
in fact be chosen equal to a, we use the sharp result of theorem 2.2. Set 
g = d@ (r s ~ 2 f R ) for some /?' < m + s — 2. Let z be an arbitrary point of B e 
and Br(z) be a ball with center 2 and radius R such that Br C -B a . Using the 
inequality (|37|), we obtain 

/ < C / ^^d/i. <C [ r^dn = C R 2+a 

Jb r (z) Jb r (z) Jb r (Q) 

(cf. [^5| p. 159 for the second estimate). Applying now Lemma 2.2, we conclude 
that g £ C a {B e ), whence r s ~ 2 f R e C m+s - 2 ' Q (£ e ). 

(iii) We have 

d P {rf R )=rdPf R + h, 
where, with certain constants Cp^pi, 

p^P'+y<p 

Note that rd' 3 / £ C a (B a ), since r is Lipschitz continuous. 
Using the bound 

\d Y r\ < Cr~^ +1 , 

and the bound ( |36|) we obtain 

< Cr a . 



Thus /i € W x ' v for all p, whence, by theorem £T], /i e C a for a < 1. 
The following function spaces will be important for us. 
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Definition 3.7 For m € No and < a < 1, we define the space E m ' a (B a ) as 
the set E m a {B a ) = {/ = fx + r/ 2 : /i,/ 2 G C m '"(S a )}. Furthermore we set 
^°°(Ba) = {/ = /i+r/ 2 : fi,h€C°°(B a )}. 

Note that the decompositions above are not unique. If / = fx + r/ 2 , /i, / 2 G 
C m ^(B a ) then also/ = / x +r/ fl + r(/ 2 - f R ) with A + rf R , f 2 - f R G C m ' a (B a ) 
if /r is given as in lemma Obviously, E co (B a ) C E m ' a (B a ) for all m G N . 
The converse is not quite immediate. 

Lemma 3.8 Iff G E m > a (B a ) for all m G N Qj i/ien / G E°°(B a ). 

Proof: Assume that / G E m ' a (B a ) for all m. Take an arbitrary m and write 
/ = f\ + rfi with /i,/ 2 G C m ' a (B a ). To obtain a unique representation, we 
write fx and / 2 as the sum of their Taylor polynomials of order m and their 
remainders, 

m m 

/i = EPi+/fl- /2 = ^P- + /I, (38) 

j=0 3=0 

with p],Pj G and = 0(r m+a ). From this we get the representation 

m rn— 1 

f = (EPi+r y Efij) + fR, (39) 

3=0 3=0 



where f R = f R + r(f 2 R +pf n ) G C m ' a (B a ) and / fl = 0(r m+a ) by lemma 3.6 



This decomposition is unique: if we had / = 0, the fast fall-off of f R at the origin 
would imply that the term in brackets, whence also each of the polynomials and 
f R , must vanish. 

Since m was arbitrary, we conclude that the function / determines a unique 
sequence of polynomials Pj, j G No as above. By Borel's theorem (cf. 
there exists a function u 2 G C°°(B a ) (not unique) such that 

m 

v 2 - J2p 2 J = 0{r m+1 ), m G N . (40) 

3=0 

We show that the function vx = f — ru 2 is C m ^ 1 (B a ) for arbitrary m, i.e. 
«i G C°°(B a ). Using @, we obtain 

(m \ / 7ii— 1 \ 

The first term is in C m ' a (B a ) by the observations above, the second term is in 
C m a (B a ) by @ and lemma |J. ■ 

While we cannot directly apply elliptic regularity results to these spaces, they 
are nevertheless appropriate for our purposes. This follows from the following 
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observation, which will be extended to more general elliptic equations and more 
general smoothness assumptions in theorem 3.11 and in appendix^. 
If u is a solution to the Poisson equation 

Am = — . 

r 

with / e E ma (B a ), m > 1, then u E E m+1 - a (B a ). This can be seen as follows. 
If we write J = fi+r f 2 G E m > a {B a ) in the form 

/ T m 

- = 1" JRi 

r r 

where T m is the Taylor polynomial of order m of /i , the remainder is seen 
to be in C m -^ a (B a ) by Lemma 

By Lemma |3.5| , there exist a polynomial T m such that 

A(rT m ) = ^. 

r 

Then ur = u — rT rn satisfies Aur — fn and theorem |2.6| implies that ur G 
C m+1 ' a (B a ), whence u e E m+1 ' a {B a ). 

To generalize these arguments to equations with non-constant coefficients 
and to non-linear equations we note the following observation. 

Lemma 3.9 For f,g<E E m ' a (B a ) we have 

(i) f + g e E m - a (B a ) 

(ii) fg&E m ' a (B a ) 

(iii) Iff^Oin B a , then 1/f e E m ' a {B a ). 
Analogous results hold for functions in E°°(B a ). 

Proof: The first two assertions are obvious, for (iii) we need only consider a 
small ball B e centered at the origin because r is smooth elsewhere. If / = 
fi + rf2, A,/2 G C m - a (B a ), we have 1/f = Vl + rv 2 with 

/i -h 



(/l) 2 -rW - (/l) 2_ r 2 (/2)2 - 

These functions are in C m,a (B t ) for sufficiently small e > because our as- 
sumptions imply that /i(0) ^ 0. The E°°(B a ) case is similar. ■ 

We consider now a general linear elliptic differential operator L of second 
order 

L = ct'didj + b l d, + c. (42) 
It will be assumed in this section that 

a ij , b\ ceC°°(B a ). (43) 



18 



We express the operator in normal geodesic coordinates x % with respect to 
a u , centered at the origin of B a , such that 

with 

fi« = 0(r 2 ), (44) 

and 

^■a y = 0, i£B„. (45) 
For the differential operator L, given by 

Lu = a v (x)didjU + b l {x)d i u + c(x)u, 

we find 

Lemma 3.10 Suppose p G P m . TTien t/ie function U defined by L(r s p) — 
r s ~ 2 U is C°° and satisfies U = 0(r m+1 ). If in addition b z = 0{r) (as in the 
case of the Yamabe operator L^), then U — 0{r m+2 ). 

Proof: A direct calculation, observing (|45|), gives 

U = a ij (s6ijp + r 2 djdip) + b^sxip + r 2 dip) + cr 2 p 7 

which entails our result. ■ 

In the following we shall use the splitting L = A + L, where A is the flat 
Laplacian in the normal coordinates x l . 

Theorem 3.11 Let u £ Wi^(B a ) be a solution of 

Lu = r s - 2 f, 

where L is given by with ffify and s € Z, p > 1. 

(i) Assume s = 1 and f G E m ' a (B a ). Then u £ E l ' a ' {B a ), < a 1 < a, if 
m = andue E m+1 > a (B a ) if m > 1. Iff e E°°{B a ), then u G E°°(B a ). 

(ii) lfs<0, \s\ is odd, f G C m > a {B a ) with m + s - 1 > 0, and f = 0{r s °), 
with sq + s > 0; then u has the form 

m 

u = r s ^2 u k + Ur, 

k=s 

where u k G V m , Lu R = 0{r m+a+s - 2 ), and u R £ C 1 ' 01 ' '(B a ), < a' < a, if 
m + s - 1 = 0, u R G C m+s > a (B a ) if m + s - 1 > 1. 

If f G C°°{B a ), then then u can be written in the form u — r s v\ + v 2 with 
v 1 ,v 2 G C°°(B a ), «]. = 0(r s °), and 

L{r s Vl )=r s - 2 f + 9 R , L(v 2 ) = -9 R , 

where 9 R G C°°{B a ) and all its derivatives vanish at the origin. 
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Proof: In both cases we write / = T m + fn with a polynomial T m of order to, 

m 

Case (i): m = and f — fi + r/2, where /i,/2 are in C m ' a (B a ). We define 
T m to be the Taylor polynomial of order to of /x- 

Case (ii): mo = so and we define T m to be the Taylor polynomial of order m of 
/• 

We show that u has in both cases the form 



u = r s ^ u k + UR, 

k—rrtQ 



(46) 



with Uk € Vk and ur S C m+S,a (i? a ). For this purpose lemma 3.5 will be used 
to determine the polynomials Uk in terms of tk by a recurrence relation. The 
differentiability of ur follows then from lemma 3.6, elliptic regularity, and the 
elliptic equation satisfied by ur. 

The recurrence relation is defined by 



A(r s u mo ) = r s -H mo , A(r s u k ) =r s - 2 (t k -U[ h) ), k = m Q + l,...,m 



(47) 



where, given u rno , . . . ,Uk-i, we define uj® as follows. By Lemma ( |3.10|) the 
functions 



fe-i 

uW =r -<+*L(r' 
which will be defined successively for k — to + 1, 



l/W _ 0(r m ° +1 ). Thus we can write by lemma (A.l) 



(48) 

, to + 1, are C°° and 



where — 0(r m+a ) and U^' G Vj denotes the homogenous polynomial of 
order j in the Taylor expansion of . 

By Lemma (3.5) the recurrence relation ( fl7| ) is well defined for the cases (i) 
e that 

U\ k ' ] = u\ k) , to + 1 < j < k if k < k' < to+ 1, (49) 



j=m + l 



(fe) 



and (ii). Note that 



because we have by Lemma 3.10 



fc'-i 



Xj(k') _ r/W = r -s+2 L{jS Uj ) = 0(r k+1 ). 

j=k 
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With the definitions above and the identity (|4^), which allows us to replace 
^fc by j the original equation for u implies for the function ur defined 



by (46) the equation 



Lun = r s 



U 



(m+l) 



fa 



(50) 



Case (i): We use lemma (3.6), (i) to conclude that the right hand side of 
equatio n (pp| ) is in L p (B e ) if m = and in C m ~ 1 ' a (B e ) if m > 1. Now theorems 
2% and |2.6| imply that u R G C 1 '"' {B a ), a' < a, if m = and u fl G C m+1 ' Q (B Q ) 
if m > 1. For the i? 00 case we use lemma 3.8. 

Case (ii): By our procedure we have 

Lur = 0(r m+a+s - 2 ). 

If 7?i + s — 1 = we use part (i) of lemma |3.6| to conclude that the right 
hand side of equation ( p0| ) is in L p (B a ) and theorems 2.5, 2.1 
Ui? G C 1,a (B a ). If m + s — 1 > 1 we use part (ii) of lemma 3.6 
the right-hand side of equation ( |50| ) is in C m+s ~ 2:01 (B a ). Elliptic regularity then 
implies that ur £ C m+s, a (B a ). The C°° case follows by an analogous argument 
as in the proof of lemma 3^, since the polynomials u mo , . . . , u m obtained for an 
integer m! with m! > m coincide with those obtained for m, i.e. the procedure 
provides a unique sequence of polynomials Uk, k = mo, ... M 



to conclude that 
to conclude that 



More general expansions, which include logarithmic terms, have been studied 
(i n a so mewhat different setting) in ]3C[|, where results similar to those given 
in 3.11 have been derived. Definition 3.7 is taylored to the case in which no 



logarithmic terms appear and leads to a considerable simplification of the proofs 
as well as to a more concise statement of the results as compared with those 
given in ]30[. 



Corollary 3.12 Assume that the hypotheses on u of theorem 3.11 are satisfied. 
Let 8q be a distributional solution of L6q = —AirSi in B a . Then we can write 



% = r 1 u 1 + u 2 , 



(51) 



withui,u 2 e C^iBa), ui(0) = 1, Lir- 1 ^) = -4nSi + 6 R> where 9 R G C* 00 ^) 
and all its derivatives vanish at i. In the particular case of the Yamabe operator 
Lh with respect to a smooth metric h we have u\ = 1 + 0(r 2 ). 



Proof: Using that A(r 1 ) — —Audi in B a , we obtain for u 

Lu = -Ur- 1 ). 



l/r 



(52) 



By lemma 3.10 we have L^r' 1 ) = r- A U, with U £ C°°(B a ) and U = 0(r). Our 



assertion now follows from theorem f3.ll . For the last assertion we use that in 
the case of Lh we have U — 0(r 2 ). ■ 



We note that the functions u%, ui are in fact analytic and 6 R = if the 
coefficients of L are analytic in B a (cf. p4j). 
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3.3 Proof of Theorem 



Ol 



There exists a unique solution 9 = 9 + u of the equation (g) with 9 as lemma 
3.2 and u as in theorem 3.4. Since the operator L b satishes the hypothesis of 
corollary 3.12 we can write on B a 

O = — +w, (53) 
r 

where u\, w are smooth functions and L) l (r~ 1 ui) = 9r on B a \ i, with 9r as 
described in corollary |3.12[ 

Given the solution u = u(x), we can read equation (|32|) in B a as an equation 
for u 



L h u 



/(*) 



(54) 



with / considered as a given function of x 

r 8 ^ ab ^ ab 



/(*) 



8(r9 



ru 



By hypothesis r 8 ^ ab ^ ab G E°°{B a ), by corollary |1| we have r 9 a € S°°(S a ), 
and by theorem |1| the solution u is in C a { B a ) C E"- a (B a ). By lemma |^9| 
wc thus have / G E°' a (B a ). By theorem 3.11 equation (|54|) implies that 
u S S 1 '"'^), < a' < a which implies in turn that / G E 1 ' 01 ' {B a ). Re- 
peating the argument, we show inductively that u, whence / is in E m,a (B a ) 
for all m > 0. Lemma 3.8 now implies that u G E°°(B a ). ■ 



3.4 Solution of the Hamiltonian Constraint with Logarith- 
mic terms 

The example 

A(logWi m ) = r- 2 h m (2ni + 1), (55) 

h m G H-m, shows that logarithmic terms can occur in solutions to the Poisson 
equation even if the source has only terms of the form r s p with p G V m ■ This 
happens in the cases where the Laplacian does not define a bijection between 
r s+2 V m and r s V m , cases which are excluded in lemma |3.5| . We shall use this 
to show that logarithmic terms can occur in the solution to the Lichnerowicz 
equation if the condition r 8 "^> ab ' l S ab G E°°(B a ) is not satisfied. Our example will 
be concerned with initial data with non-vanishing linear momentum. 

We assume that in a small ball B a centered at i the tensor \l/ ab has the form 

where ^p is given in normal coordinates by ( |76| ) and = 0(r~ 3 ) is a tensor 
field such that ^ ab satisfies equation (0). The existence of such tensors, which 
satisfy also 

8o 7 *ab* ab eL«(S), q> 2, (57) 
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and, by lemma 



4.13 



ah 



tp + ril> R 



B„ 



(58) 



will be shown in section 4.3. Here the function tp R is in C a (B a ) and ip is given 
cxplicity by 



tp = r^pyf/ = c + r z h 2 



(59) 



with c = ±| P 2 , P 2 = pip i} and 



h 2 = -r 2 (3(P l ni 

o 



p2 



(60) 



where, in accordance with the calculations in section 4.3, the latin indices in 



the expressions above are moved with the flat metric. We note that h 2 G H 2 
and if) is not in continuous. The tensor ^ ab satisfies condition (|^) and the three 
constants P l given by (76) represent the momentum of the initial data. Since 
VPp 6 is trace-free and divergence-free with respect to the flat metric, we could, 
of course, choose h a f, to be the flat metric and ^> R h — 0. This would provide one 
of the conformally flat initial data sets discussed in Jll| . We are interested in a 
more general situation. 



Lemma 3.13 Let h a {, be a smooth metric, and let \E ,ab be given by ft5q ) such that 
conditions §5\) and $5f\ ) hold. Then, there exists a unique, positive, solution to 
the Hamiltonian constraint |s|). In B a it has the form 



- + l m + ^ r ((9(PX) 2 -33P 2 ) 
7 



Wi 

r 

+ m mr ' 2 [jP 2 + 1^{p 1 ^) 2 -p 2 )^^+u r , (6i) 

where the constant m is the total mass of the initial data, w\ is a smooth function 
with wi = 1 + 0(r 2 ), and u R S C 2 ' a (B a ) with u R (0) = 0. 

Since w\ is smooth and u R is in C 2 ' a (B a ), there cannot occur a cancellation 
of logarithmic terms. For non-trivial data, for which m ^ 0, the logarithmic 
term will always appear. In the case where /i a j, is flat and = an expansion 
similar to ( |si|) has been calculated in |2^| . 

Pro of: The existence and uniqueness of the solution has been shown in section 
3.1 . To derive ( |6l"| ) we shall try to calculate each term o f the expansion and to 
control the remainder as we did in the proof of theorem 3.11. However, lemma 
3.5 will not suffice here, we will have to use equation (|55). 
By Corollary 3.12 we have 



r 
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with w,wi £ C°°(B a ) and W\ = 1 + 0(r 2 ). By theorem 3.4 the unique solution 



u of equation ( p4| ) is in C a (B a ). Equation ( |54| ) has the form 

ib 1 

LhU= (- +4> R ) f(x,u) with f(x,u) = -- ■ — - — ■ rry. 

r (wi + r(u + w))' 

By m = 2 (u(0) + w(0)) is given the mass of the initial data. Since u £ C a (B a ), 
we find 

f = -l + \mr + f R , f R = 0(r 1+a ). (62) 



If we set 



we find from @, @ that 



t/» 7 

A«i = -(-l+-mr), (64) 



r v 2 



and that i> = u — iti satisfies 



L h v = ib R (-l + lmr)-L hUl + ^. (65) 
2 r 

We shall show that this equation implies that the function v is in C 2,a {S). 
Since i> R £ C a {B a ), the first term on the right hand side of © is in C a (B a ). 
By a direct calculation (using that the coefficient b % of Lh is O(r)) we find that 
LhU\ £ W 2,p (B a ), p < 3, which implies by the Sobolev imbedding theorem 
that L^ui £ C a (B a ). The third term on the right-hand side of (pq ) is more 
delicate, because it depends on the solution v. From theorem [3.4| and equation 
© we find that v £ C a (B a )f]W 2 ' p (B a ), where, due to @, we need to assume 
p < 3. However, since /r = Q(r 1+a ) and t/» is bounded, the function ipfn/r is 



bounded and the right hand side of ( |65| ) is in L°°(B a ). Thus theorem 2.5 implies 
that u € W 2 ' p (B a ) for p > 3, whence u € C 1 ' a (B a ) by the Sobolev imbedding 
theorem. It follows that we can differentiate fn, considered as a function of x, 
to find that d t f R = 0{r). It follows that is in W^ v {B a ), p > 3, and 

thus in C a (B a ). Since the right hand side of ( |65| ) is in C a (B a ), it follows that 



4 The Momentum Constraint 

4.1 The Momentum Constraint on Euclidean Space 

In the following we shall give an explicit constructing of the smooth solutions to 
the equation d a ^ ab = on 3-dimensional Euclidean space E 3 or open subsets of 
it. Another method to obtain such solutions has been described in 0, multipole 
expansions of such tensors have been studied in [H . 
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Let i be a point of E 3 and x k a Cartesian coordinate system with origin i 
such that in these coordinates the metric of E 3 , denoted by S ab , is given by the 
standard form Ski- We denote by n a the vector field on E 3 \ {i} which is given 
in these coordinates by x k /\x\. 

Denote by m a and its complex conjugate fh a complex vector fields, defined 
on E 3 outside a lower dimensional subset and independent of r = \x\, such that 

m a m a = m a rh a = n a m a = n a m a = 0, m a fh a — 1. (66) 

There remains the freedom to perform rotations m a — > e l *m a with functions / 
which are independent of r. 

The metric has the expansion 

Sab = n a n b + fh a m b + m a fh bl 

while an arbitrary symmetric, trace-free tensor ty ab can be expended in the form 

r 3 *ab = £(3n a n & - 5 ab ) + \f2^n (a m h) + V2fjin^ a m b ) + 

\iim a m b + fj, 2 m a fhb, (67) 

with 

£ = ^r 3 ^> ab n a n b , i h = V2r 3 f ab n a m b , fi 2 = r 3 ^ ab m a m b . 

Since ^> a b is real, the function £ is real while rji, \ii are complex functions of 
spin weight 1 and 2 respectively. 
Using in the equation 

d a V ab = 0, (68) 

the expansion (|67]) and contracting suitably with n a and m°, we obtain the 
following representation of ( pH ) 

4rd r £ + Br)! + dfji = 0, (69) 



rd rm + 8 M2 - S£ = 0. (70) 

Here d r denotes the radial derivative and S the edth operator of the unit two- 
sphere (cf. ]32| for definition and properties). By our assumptions the differen- 
tial operator S commutes with d r . 

Let s Yi m denote the spin weighted spherical harmonics, which coincide with 
the standard spherical harmonics Yi m for s = 0. The s Yim are eigenfunctions of 
the operator S3 for each spin weight s. More generally, we have 

[s lm) { l) (l-8-p)\ (1 + S)\ S lm - { ' 



25 



If n a denotes a smooth function on the two-sphere of integral spin weight 
s, there exists a function /i of spin weight zero such that i] s = d s n. We set 
n R = Re(ry) and rf — ilm.(rj), such that n = r] R + i] 1 , and define 

such that rj s = n R + n I s . It holds 



& V R = Wri R , SV = -g s 7? J . 

Using these decompositions now for 771 and /^2, we obtain equation ( |69|) in 
the form 

2rcU = -BBr) R . (72) 

Applying 9 to both sides of equation ( |70| ) and decomposing into real and imag- 
inary part yields 

rdrddr] 1 = -6 2 9V (73) 
2rd r (rd r + S3£ = 8 2 g 2 /A (74) 



Since the right hand side of ( |72[ ) has an expansion in spherical harmonics 
with I > 1 and the right hand sides of ([73]), (|74|) have expansions with I > 2, 
we can determine the expansion coefficients of the unknowns for / = 0, 1. They 
can be given in the form 



with 



i = A + rQ + -P, i] 1 = iJ + const., r\ R = rQ --P + const, 
r r 



P=^P a n a , Q=^Q a n a , J = 3J a n a , (75) 



where A,C a ,Q a , J a are arbitrary constants. Using d67j), we obtain the corre- 
sponding tensors in the form (cf. (JTl|) 



_ _3_ ^_pa n 6 _ p b n a _ (gab _ 5n a n fc^ p c n j j ^g) 
*j ^ 4( nOe60<Ij ° n * + n h e acd J c n d ), (77) 



^> = -3 (3n a n b - <5 ob ), (78) 

*q = ^2 (Q"" 6 + Q bn<1 - ( sab - nClnb ) Q Cjl -) ■ ( 79 ) 

We assume now that £ and r\ have expansions in terms of in spherical 
harmonics with I > 2. Then there exists a smooth function A2 of spin weight 2 
such that 
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Using these expressions in equations ( [72] ) - ( f74| ) and observing that for smooth 
spin weighted functions fj, s with s > we can have dfi s = only if ji s = 0, we 
obtain 

5?7 fl = -2r<9 r 8A^, g 2 // = -^ r A^, 

o 52 ^ = 2rd r {rd r \f) - 2A^ + 33A? . 

We are thus in a position to describe the general form of the coefficients in the 
expression ( |67| ) 

£ = d 2 \X + A + rQ + -P, (80) 
r 

r/i = -2rd r B\$ + B\i + rBQ- -BP + iBJ, (81) 

r 

/i 2 = 2rd r (r d r \%) ~ 2 \* + 9§Af - r d r \[. (82) 

Theorem 4.1 Let A oe an arbitrary smooth, complex, function in B a \{i} C IE 3 
with < a < oo, and set A2 = S 2 A. T/ien i/ie tensor 

yab _ ^ab + ^ah + ^ah + ^ab + ^cri (gg) 

satisfies D ai i ^ — Q in B a \ {i}. Here the first four terms on the right hand side 
are given by (f^jj - f|76]j while is is obtained by using in only the part 
of the coefficients j&uT - ftsdj ) which depends on A2. Conversely, any smooth 



solution in B a \ {i} of the equation above is of the form (85) 



Obviously, the smoothness requirement on A can be relaxed since "f^ 6 G 
C 1 (-B a \ {i}) if A S C 5 (B a \ {i}). Notice, that no fall-off behaviour has been 
imposed on A at i and that it can show all kinds of bad behaviour as r — > 0. 

Since we are free to choose a, we also obtain an expression for the general 
smooth solution on E 3 \ {i}. By suitable choices of A we can construct solutions 
vp^ b which are smooth on E 3 or which are smooth with compact support. 

Given a subset S of M 3 which is compact with boundary, we can use the 
representation ( |S3|) to construct hyperboloidal initial data (^l)) on S with a 
metric h which is Euclidean on all of S or on a subset U of S. In the latter case 
we would require 'J'^ 6 to vanish on S \ U. In the case where the trace-free part 
of the second fundamental form implied by h on dS vanishes and the support 
of ^> ab has empty intersection with dS the smoothness of the corresponding 
hyperboloidal initial data near the boundary follows from the discussion in ( § ) . 
Appropriate requirements on h and *& ab near dS which ensure the smoothness 
of the hyperboloidal data under more general assumptions can be found in ( Q ) . 

There exists a 10-dimensional space of conformal Killing vector fields on E 3 . 
In the cartesian coordinates x % a generic conformal Killing vector field £q has 
components 

C a = k j (2 Xj x l - Sj l xix l )+ e l jk S j x k + ax 1 + q l , (84) 
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where k l , S 11 , q l are arbitrary constant vectors and a an arbitrary number. 
In terms of the "physical coordinates" y % — x l /\x\ 2 , with respect to which i 
represents infinity, we see that k l , S l , q l , and a generate translations, rotations, 
"special conformal transformations" , and dilatations respectively. 

For < e < a we set S e = {\x\ = e] and denote by dS e the surface element 
on it. For the tensor field ^ ab of (|33|) we obtain 



^- / * ab n a ^ b dS e = (P a k a + J a S a + Aa + Q a q a ). 



(85) 



We note that the integral is independent of e and, more importantly, indepen- 
dent of the choice of A. Thus the function A does neither contribute to the 
momentum 

P a = -!-lim/ r 2 y bc n b (2n c n a -S ca )dS t , (86) 
nor to the angular momentum 

J a = -!- lim / r V bc n b e cad n d dS e , (87) 

of the data. 

If we use the coordinates x % to identify E 3 with R 3 and map the unit sphere 
S 3 C K 4 by stereographic projection through the south pole onto K 3 , the point 
i, i.e. x l = 0, will correspond to the north pole, which will be denoted in 
the following again by i. The south pole, denoted in the following by i' , will 
represent infinity with respect to the coordinates x % and the origin with respect 
to the coordinates y 1 . We use x 1 and y l as coordinates on S 3 \ {%'} and S 3 \ {i} 
respectively. If h° is the standard metric on S 3 , we have in the coordinates x l 

I , „2\ 1/2 



h° kl = e- 2 s kl , e=^^-j . (88) 

We assume that the function A is smooth in E 3 \ {i} and set ^ ab = 6 10 * ab 



with f ah as in (B§. Then, by (111), <M 



D a * ab -0 in S 3 \{i,i'}, (89) 

where D a denotes the connexion corresponding to h°. The smoothness of ty ab 
near i can be read off from equations ( |67| ) and (|8(]) - (|2|). In order to study 
its smoothness near i' we perform the inversion to obtain the tensor in the 
coordinates y l . It turns out that we obtain the same expressions as before if we 
make the replacements 

n % — > -n\ m % — > m\ 
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P a ^Q a , J a ^-J a , A^A, Q a ^P a . 



Thus the tensors ( |76| ) - (pg), the first two of which are the only ones which 
contribute to the momentum and angular momentum, are singular in i' as well 



as in i. Observing again the conformal covariance (113), we obtain the following 
result. 

Corollary 4.2 The general smooth solution of the equation D a ^ ab = on 
S 3 \{i,i'} with respect to the metric ui 4 h° , where h° denotes the standard metric 
on the unit 3-sphere and uj S C°°(S 3 ), uj > 0, is given by * ofc = (uj^ 1 9) w * ab 
with ^/ a (, as in ( f&| j and A € C°°(E 3 \{i}). If we require the solution to be bounded 
near i' ( in particular, if we construct solutions with only one asymptotically flat 
end), the quantities P a , J a , A, Q a must vanish. 

We can now provide tensor fields which satisfy condition (Thy) and thus prove 



a special case of theorem 1.2 



Theorem 4.3 Denote by * ab a tensor field of the type If rX e E co {B a ) 

and P a = 0, then r s ^ ab ^ ab £ E°°(B a ). 

Only the part of A which in an expansion in terms of spherical harmonics 
is of order I > 2 contributes to ^ a b- We note that the condition rA G E°°(B a ) 
entails that this part is of order r. The singular parts of are of the form 

©-©• 

For the proof we need certain properties of the S operator. If m a is suitably 
adapted to standard spherical coordinates, the S operator, acting on functions 
of spin weight s, acquires the form 



5r/ s = — (sin( 



d id 



89 sin ( 



((sin0)-%). (90) 



The operator 33 acting on functions with s — is the Laplace operator on the 
unit sphere, i.e. we have the identity 

33/ = r 2 Af - y.r'i),,),f - 2x i d i f, (91) 

where A denotes the Laplacian on R 3 . The commutator is given by 

(33 - 33)?7 S = 2sn s . (92) 

From this formula we obtain for q S N by induction the relations 

(§a* _ m)n s = {2qs + (q- l)q)W- l ri s , (93) 

(33 9 - B q B)rj s = (-2qs + (q- l)q)3 9 - V (94) 

In particular, it follows that the functions rj and fj,, which satisfy 3?7 = r\\ and 
3 2 /i = /i2, are given by 

7} = -2rd r (33A fl + 2\ R ) + 33A' + 2X 1 + rQ - — + U, (95) 

r 

H = 2rd r {rd r \ R ) + m\ R - 2rd r X I . (96) 
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Lemma 4.4 Let /, g G E°° be two complex functions (spin weight zero). Then, 
the functions d q fd q g, g£N, of spin weight zero are also in E°° . 

Remarkably, the statement is wrong if we replace E°° by C°° : the calculation 
in polar coordinates, using (p0|), gives 

Qx 1 Bx 3 = -x 1 x 3 + irx 2 . 

Proof: For q = 1 the proof follows from two identities. The first one is a simple 
consequence of the Leibnitz rule 

dfBg + 8/g 5 = m(fg) - fm g - g mf. (97) 

Since §3 maps by ( |9l| ) smooth functions into smooth functions, it follows that 
(dfBg + Bfdg) G E°° if /, g G E°° (here we can replace E°° by C°°). 
The other identity reads 

9/8.9 - 57% = 2 i r e^ k x l djf 8 k g. (98) 

It is obtained by expressing ( |90| ) in the Cartesian coordinates x l . Important for 
us is the appearance of the factor r. A particular case of this relation has been 
derived in [ p>7j. It follows that (BfBg - Bfdg) G E°° if /, g G E°° . Taking the 
difference of (|97|) and ( |98| ) gives the desired result. 

To obtain the result for arbitrary q, we proceed by induction. The Leibniz 
rule gives 

B q+1 fB q+1 g = Bd{d q fB q g) - m q fdd q g - B q fdd q+1 g - d q gdm q f. 

The induction hypothesis for q and ( pl| ) imply that the first term on the right 
hand side is in . The factors appearing in the following terms can be written 



by (|93j) and (|9j) in the form 

m q f = d q -\mf + q(q - 1)/), m q g = Bi-^mg + q(q - l)g), 

m q+1 g = B q (mg + q(q - l)g), md q f = W(BBf + q(q - 1)/). 

Since the functions in parenthesis are, by (|9l|), in E°° , the induction hypothesis 
implies that each of the products is in E°° . ■ 



Proof of theorem 4.3: In terms of the coefficients (pOf) - (fe2) we have 



r s y ab V ab = r 2 (2 M2 /2 2 + 2rnfk + 3^) ■ (99) 
Since rX is in E°°(B a ), equations (Ba) and (|9^) imply that r^,rn,r^ G E°°. 



The conclusion now follows from lemma 4.4 
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4.2 The general case: existence 

The existence of solution to the momentum constraint for asymptotically flat 
initial data has been proved in weighted Sobolev spaces (cf . O] , [jl5| , Q , 
and the reference given there) and in weighted Holder spaces |13| . The existence 
of initial data with non trivial momentum and angular momentum and the role 
of conformal symmetries have been analysed in some detail in j^] . In this section 
we will prove existence of solutions to the momentum constraint with non-trivial 
momentum and angular momentum following the approach of Q . We generalize 
some of the results shown in |J to metrics in the class (^lj) . The results of the 
previous section will be important for the analysis of the general case. 

We use the York splitting methods to reduce the problem of solving the 
momentum constraint to solving a linear elliptic system of equations. Let the 
conformal metric h on the initial hypersurface S be given. We use it to define 
the ovcrdctcrmincd elliptic conformal Killing operator Ch, which maps vector 
fields v a onto symmetric ft,-trace-free tensor fields according to 

(£ h v) ab = D a v b + D b v a - ~ h ab D c v c , (100) 

o 

and the underdetermined elliptic divergence operator Sh which maps symmetric 
/i-trace-free tensor fields $ ab onto vector fields according to 

(5 h $) a = D b <P ba . (101) 

Let $ afc be a symmetric /i-trace-free tensor field and set 

V ab = <$> ab -(£ h v) ab . (102) 

Then ^ ab will satisfy the equation D a ^f ab = if the vector field v a satisfies 

L h v a = D b <S> ab , (103) 

where the operator is given by 

L h v a = D b (C h v) ab = D b D b v a + ^D a D b v b + R a b v b . (104) 

Since ^ ab (C h v) ab = 2D a (^ ab v b ) - 2 (5 h V) a v a for arbitrary vector fields v a and 
symmetric /i-trace-free tensor fields ^ ab , Ch has formal adjoint L\ = —2Sh- 
Thus 

L h = ~C* h ojC h , (105) 

and the operator is seen to be elliptic. 

Provided the given data are sufficiently smooth, we can use theorem 2.7 to 



show the existence of solutions to (103) 



Lemma 4.5 (Regular case) Assume that the metric satisfies pj^ j and <& ab is 
an h-trace-free symmetric tensor field in W 1 ' P (S), p > 1. Then there exists a 
unique vector field v a € W 2 > P (S) such that the tensor field * afc = $ ab - (C h v) ab 
satisfies the equation D a i 5> ab = in S. 
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Proof: The well known argument that the condition of theorem 2.7 will be 
satisfied extends to our case. Assume that the vector field £ a is in the kernel of 
Lh, i-e. 

L h C = in S. (106) 
Since the metric satisfies (PH), elliptic regularity gives 

C e C 2 ' a (S). (107) 



This smoothness suffices to conclude from ( |105| ), fllOq ) that = — 2 (£, L^)^ = 
(Ch£„ChCjL 2 , whence 

(ChOab = 0. (108) 

This implies for an arbitrary symmetric, /i-trace-free tensor field $ ab £ W 1,p (5), 
p > 1, the relation = (Ch£, &)l 2 = — 2 (£, <fyi$)i 2 5 which shows that the Fred- 
holm condition will be satisfied for any choice of $ afc in (103). ■ 

We call the case above the "regular case" because the solution still satisfies 
the condition fy ab £ W X ' P (S). While this allows us to have solutions diverging 
like 0(r~ 1 ) at given points, it excludes solutions with non- vanishing momentum 
or angular momentum. 

We note that by fllOSD the kernel of L/j consists of conformal Killing fields. 
Let £ a be such a vector field. Using (107) and lemma |A.1| we find that we can 



write in normal coordinates centered at the point i of S 

e fc =£o +0(r 2+a ), (109) 
where £q is the "flat" conformal Killing field (Q) with coefficients given by 

k a =\D a D h ?(i), S a = e a bc D b H c (i), <f = e(i), a=^D a C(t). (HO) 
6 3 

Since S is connected, the integrability conditions for conformal Killing fields (cf. 
[pi]) entail that these ten "conformal Killing data at i" determine the field £ a 
uniquely on 5*. 

With a conformal rcscaling of the metric with a smooth, positive, conformal 
factor lo 

h a h — * h' ab = uj 4 h ab , (111) 

which implies a transition of the connection D a — > D' a , we associate the rescal- 
ings 

y ab _> ^' ab = cj- 10 * Qb , £ /a = £°, (H2) 

for ft-trace free, symmetric tensor fields ^ ab and Killing fields £ a . Then the 
conformal Killing operator and the divergence operator satisfy 

(C h ,v) ab = cj-\£ h v) ab , D' a * ,ab = uj- 10 D a V ab . (113) 
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If we write u> = e^, the conformal Killing data transform as 

k' a = k a + 2aD a f(i) + e abc S b D c f(i) + q c D a D c f{i), (114) 

S' a = S a + 2e abc q b DJ(i), (115) 

a' = a + q a D a f(i) (116) 

q' a = q a . (117) 

A vector field £ on S will be called an "asymptotically conformal Killing 
field" at the point i if it satisfies in normal coordinates x k centered at i equa- 



tion (109) with (|84|). While any conformal Killing field is an asymptotically 
conformal Killing field, the converse need not be true. Let £ a be an asymptoti- 
cally conformal Killing field and consider the integral 

/? = -Lfim / V ab n a t; b dS e , (118) 

where in the coordinates x k the unit normal n to the sphere dB e around i 
approaches x fc /|a;| for small e. As shown in the previous section, the integral 
vanishes if the tensor *5> ab is of order r _1 at i. If, however, ^ ab is of order r~ n , 
n — 2,3,4, at i, the integral gives non-trivial results and can be understood in 
particular as a linear form on the momentum and angular momentum of the 
data. 

We recall two important properties of the integral 1^. The first is the fact 
that the integral 1^ is invariant under the rescalings (111), (112). The second 
property is concerned with the presence of conformal symmetries. Let $° be 
an arbitrary /i-trace-free tensor field and v a an arbitrary vector field. Using 
Gauss' theorem, we obtain 

2/ 1> ab n a v b dS e = -2 f v a D b $ ab dn- [ (C h v) ab ^ ab dfi, (119) 

JdB e JS-B e JS-B e 



with the orientation of n k as above. If we apply this equation to (118), the first 
term on the right hand side will vanish if Da^ ab = 0inS'\{i}, while the second 
term on the right hand side need not vanish if £ a is only an asymptotically 
conformal Killing vector field. However, if £ is a conformal Killing field, we get 
Jf = 0. Thus the presence of Killing fields entails restrictions on the values 
allowed for the momentum and angular momentum of the data. For this reason 
the presence of conformal symmetries complicates the existence proof. Note 
that we are dealing with vector fields £ a which satisfy the conformal Killing 



equation (108) everywhere in S; in general, a small local perturbation of the 



metric will destroy this conformal symmetry. 



Observing the conformal covariance (113) of the divergence equation, we 
perform a rescaling of the form (|25|), ([27 ) such that the metric h' has in h'- 
normal coordinates x k in B a centered at i local expression h' kl with 

h' kl - Su = 0(r 3 ), d 3 ti kl =0{r 2 ). (120) 
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In these coordinates let ^^ oi € C°°(B a \ {i}) be a trace free symmetric and 
divergence free tensor with respect to the fiat metric Ski, 

6 ik *f lat = 0, d^f lat = in B a \ (121) 

with 

*}L = 0(r- 4 ), djVfu* = 0(r- 5 ) as r - 0. (122) 



All t hese tensors have been described in theorem ( [f .l[) . Note that the conditions 
(122) are essentially conditions on the function A which characterizes the part 
3f Q. Denote by <5>f lng G C°°(S \ {i}) the /i'-trace free tensor which is 
given on B a \ {i} by 

= x mlt - J ^ h ia \ (123) 

and vanishes elsewhere. Here x denotes a smooth function of compact support 
in B a equal to 1 on B a / 2 - By our assumptions we have then 

D'A q = 0(r- 2 ) as -> 0. (124) 



Theorem 4.6 (Singular case) Assume that the metric h satisfies n2i), ujq 



denotes the conformal factor h' — loq h, $" 4 „ g is the tensor field defined 

above, and Q^eg * s a symmetric h-trace free tensor field in W 1,P (S), p > 1. 

i) If the metric h admits no conformal Killing fields on S , then there exists a 
unique vector field v a € W 2,q (S), with q = p if p < 3/2 and 1 < q < 3/2 if 
P > 3/2, such that the tensor field 

* ab = ^ + <ff eg + (£ h ,v) ab ), (125) 

satisfies the equation D a ^ ab = in S \ {i}. 

ii) If the metric h admits conformal Killing fields £ a on S , a vector field v a as 
specified above exists if and only if the constants k a , S a , a, q a (partly) charac- 
terizing the tensor field $°f„ g (cf. (8S), ( 76)-^7^)), satisfy the equation 

P a k a + J a S a +Aa+{P c L c a (i) + Q a )q a = 0, (126) 

for any conformal Killing field t; a of h, where the constants k a , S a , a, q a char- 
acterizing £ a are given by (11C). 

In both cases the momentum and angular momentum (cf. $8d[), of^ ab 
agree with those of the tensor &s b ng - These quantities can thus be prescribed 
freely in case (i). 
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Proof: Because of (124) wc can consider D' a (^ b ng + ^eg) as a function in 
L P (S), 1 < p < 3/2. In case (i) the kernel of the operator Ch> appearing in the 
equation 



D' ($ afc 

a V sin 6 



$? eg + {L h ,vY b ) = 0, 



is trivial and we can apply theorem 2.7 to show that the equation above deter- 
mines a unique vector field v a with the properties specified above. After the 
rescaling we have D a 'i> ab = by (113). 

In case (ii) the kernel of Ch> is generated by the conformal Killing fields 
= £ a of h' . If we express ( |119| ) in terms of the metric h! , the tensor field 
+ $"e„, and the vector field £' a , take the limit e — > and use equation 
5|), we find that the Fredholm condition of theorem 2.7 is satisfied if and only 
if for every conformal Killing field £ /Q of h! we have 

P a k' a + J a S' a + Aa' + Q a q' a = 0, 



$> ab 

sin g 



where the constants k u 
terms of £ /a and h'. By Jill) 



a', q' n a re given by equation (IK), exp resse d in 
(117) this condition is identical with (126). 



As shown in the previous section, we can choose $°f ns such that the corre- 
sponding momentum and angular momentum integrals take preassigned values, 
which can be chosen freely in case (i) and need to satisfy (|l26| ) in case (ii). 
These values will agree with those obtained for uu Q ~ w '$> ab due to the regularity 
properties of v a . After the rescaling the values of the momentum and the an- 
gular momentum remain unchanged because luq = 1 + 0{r 2 ). M 



We note that it is the presence of the 10-dimensional space of conformal 
Killing fields on the standard 3-sphere which led to the observation made in 



corollary (4.2). The latter generalizes as follows. 



Corollary 4.7 If S is an arbitrary compact manifold, h satisfies (piy, and we 
allow for p > 2 asymptotic ends ik, 1 < k < p, we can choose the sets of 
constants (P£ , J£, Ak,Q%) arbitrarily in the ends ik, 1 < k < p — 1. Which 



constants can be chosen at the 
admitted by h. 



id i p depends on the conformal Killing fields 



Proo f: This follows from the observation that in the case of p ends equation 
(126) generalizes to an equation of the form 



J2 Pt K + J? S l a +A t a l + (Pf L c Q ( i; ) + Qf) q l a = 0, 



i=i 



where the constants bear for given I the same meaning with respect to the point 
%i as the constants in (126) with respect to i. ■ 



The case of spaces conformal to the unit 3-sphere (S 3 , ho) is very exceptional. 
A result of Obata p3|, discussed in in the context of the constraint equations, 



35 



says that unless the manifold (S, h) is conformal to (5 3 , ho) there exists a smooth 
conformal factor such that in the rcscaled metric h' every conformal Killing field 
is in fact a Killing field. Thus the dimension of the space of conformal Killing 
fields cannot exceed 6. In fact, it has been shown in Q that in that case h' 
can admit at most four independent Killing fields and only one of them can be 
a rotation. In this situation equation (|126|), written in terms of the metric h! , 



reduces by (110) to 

J a S' a + (Q a + P h L' a b (i))q a = 0, 

since D' a ^ a = for a Killing field. The constants P a and A can be prescribed 
arbitrarily. If there does exist a rotation among the Killing fields, the equation 
above implies 

J a S' a = 0, (Q a + P b L' a b (i))q a =0. 

4.3 Asymptotic expansions near i of solutions to the mo- 
mentum constraint 



In this section we shall prove an analogue of Theorem 3.11 for the operator L/j 



defined in (104). It will be used to analyse the behaviour of the solutions to 



the momentum constraint considered in theorem (4.6) near i and to show the 
existence of a general class of solutions which satisfy condition (0). Our result 
rests on the close relation between the operator and the Laplace operator. 

We begin with a discussion on R 3 and write xi = x z , d l = di- The flat space 
analogue of L/ t on R 3 is given by 

L v k = Av k + ^d k d t v l , (127) 

where A denotes the flat space Laplacian and v k a vector field on some neigh- 
bourhood of the origin in R 3 . 

The following spaces of vector fields whose components are homogenous 
polynomials of degree m and smooth functions respectively will be important 
for us. 

Definition 4.8 Let m £ N, m > 1. We define the real vector spaces Q m , 
Qoo{B a ) by 

Q m = {v£ C°°(M 3 ,R 3 )| v l £ V m , V*Xi = r 2 v With V £ V m -i}, 



Qoo(Ba) = {v£ C°°(B a ,R 3 )\ v l x, = r 2 v with v £ C°°{B a )}. 



The following lemma, an analogue of lemma B.J: , rests on the conditions imposed 
on the vector fields above. 
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Lemma 4.9 Suppose s € Z. Then the operator Lo defines a linear, bijective 
mapping of vector spaces 



^o-r s Q m ^r s - 2 Q m , 

in the following cases: 

(i) s > 

(ii) s < 0, \s\ is odd and m + s > 0. 

Note that the assumptions on m and s imply that the vector field Lo(r s p 1 ) £ 
C°°(R 3 \{0}, M 3 ) defines a vector field in ^(K 3 , K 3 ) which represents L (r s p*) 
in the distributional sense. 

Proof: For s as above and p l € Q m there exists some p m -i S V m —\ with 

dk{r s p k ) = r s q m -i with g m _ x = sp m _ x + d k p k E V m -\- (128) 
With equation ( |34| ) it follows that 

L (rV) = r s_2 p' with 



f = s (s + 1 + 2m) p ,: + r 2 Ap J + - (s ar* + r 2 d l q m ^) E V m . 
Moreover, p 1 E Q m because p l Xi = r 2 p m -i with 

Pm-l = S (s + 1 + 2m) p m -x + Xi Ap 1 + -(s q m -l + Xi d l q m -l) € Vrn-\- 

To show that the kernel of the map is trivial, assume that ho(r s p l ) =06 
L i 1 oc (R 3 , R 3 ). Taking a (distributional) derivative we obtain 

= 5iL (rV) = ~ A(ft(rV)) = i A(r s q m -i)- (129) 



When s > or \s\ odd and m — 1 + s > we use lemma ^5| to conclude that 
<7m-i = 0. We insert this in the equation ~Lo(r s p l ) = to obtain A(r s p J ) = 
and conclude again by lemma [T^ that p 1 = 0. 

There remains the case s + m = with \s\ odd. Expanding g m _i in equation 



(129) in harmonic polynomials (cf. (pq)), we get 



= A(r- m q m ^) = ^ A(r 2fe_m /i m _i_2fe), 

0<fc<(m-l)/2 



whence, by (|34|) , 



53 r 2k ~ m - 2 (m - 2 fc) (m - 2 fc - 1) fr ro _i_ 2fe = 0. 

0<fc<(m-l)/2 
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Since this sum is direct each summand must vanish separately. Since m is odd, 
the only factor (m — 2k — 1) which vanishes occurs when 2k = m — 1, from 
which we conclude that q m -i = r" l_1 /io with a constant ho- Since equation 



(129), which reads now ho A(r *) = 0, holds in the distributional sense, it fol- 



lows that ho = 0. ■ 

Unless noted otherwise we shall assume in the following that the metric h 
is of class C°° and that it is chosen in its conformal class such that its Ricci 
tensor vanishes at i (cf. (|25|), (|27|)). By x % will always be denoted a system of 
/i-normal coordinates centered at i and all our calculations will be done in these 
coordinates. Thus we have 

h M =5 kl +0{r 3 ), d h kl =O(r 2 ). 

We write the operator in the form 

L/i = Lo + L/,, 

where, with the notation of ([22]), 

(L h v)i = h^djdkVi + ^h jk did k Vj + B> k idj v k + A> i Vj , (130) 



with 



b* i = -2 h jl r, k i - - h l f r, * ' / h k l + - ft h jk 



and A 3 j is a function of the metric coefficients and their first and second deriva- 
tives. The fields A> h B kj i are smooth and satisfy 

A^ = 0(r), B fc ^ = 0(r 2 ), (131) 

and, because x k x l Ti k j — at the point with normal coordinates x k , 

x k x i B k j i = ~r ;, h lf r^ f . (132) 
Similarly, we write the operator Ch in the form 

£h = Cq + £h- 

Lemma 4.10 Suppose p % 6 Q m . Then Lh(r s p l ) = r s ~ 2 U l with some U l <G 
Qoo{B a ) which satisfies U l = 0{r m + 3 ). 



Proof: Using ( |130| ), we calculate r s+2 hh{r s p l ) and find 

Ui = h k3 (sS kjP l + r 2 d k d jPl ) + ~h jk (sxApj + r 2 d k d lPj + 8 klPj ) 

+ B kj i (sxkPj + r 2 d kVl )+r 2 AUpj. (133) 
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Thus U l is smooth. Using fllM] ) we obtain that U t = <3(r m + 3 ), using ( p^ ) we 
find o^t/j = r 2 f with some smooth function /. ■ 

We are in a position now to prove for m = oo the analogue of point (ii) of 
theorem 3.11. 

Theorem 4.11 Assume that h is smooth, s £ Z, s < 0, |s| odd, F 1 € C°°(B a ), 
and J 1 £ Qoo(-B a ) wi£/i J 1 = O(r so ) /or some So > \s\. 
Then, ifv 1 e W?£{B a ) solves 

{Uvf = r s ~ 2 J i + F i , 

it can be written in the form 

v l =r s v\+v l 2 , (134) 
with v\ £ QociBa), v\ = 0(r s °), v\ £ C°°(B a ). 



Proof: The proof is similar to that of Theorem 3.11. For given m £ N we 
can write by our assumptions J 1 = + J R , where J R — 0(r m+1 ) and 
denotes the Taylor polynomial of J 1 of order m. Because J 1 £ Qoo(B a ), its 
Taylor polynomial can be written in the form 

m 

Tl = £ 4 with 4 G Qk- 

k=so 

We define now a function v R (depending on m) by 

m 

"'^ s ^ vi + Vr . 

k=so 

The quantities (v k ) £ Qk are determined by the recurrence relation 



L (rX 



*- 2 t* Lo(rX) = r s " 2 (4 - C/p), 



where, for given fc, the quantity U^ 1 £ Qk is obtained as follows. The function 



jj(k)i 



k-1 



-s+2 



has by lemma 4.1C an expansion 



T (k)i 



(k)i rr(k)i 
j + U R 



with U<j k)i £ Qj, 



from which we read off . By Lemma 4.9 the recurrence relation is well 
defined. 
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With these definitions, the remainder v l R satisfies the equation 



U 



(m+l)i 
R 



Jl 



F\ 



By lemma B.6 the right hand side of this equation is in C m+S 2,a (B e ). By ellip- 
tic regularity we have v R S C m+S,a (B a ). Since m was arbitrary, the con clus ion 



follows now by an argument similar to the one used in the proof of lemma ^8 



Theorem 4.11 will allow us to prove that the solutions of the momentum 
constraint obtained in section 4.2 have an expansion of the form ([l^), if we 
impose near i certain conditions on the data which can be prescribed freely. 



In definition ( |123| ) of the field $" 2 h „ g , which enters theorem ( |4.6| ), we assume 
t ^££ ot is of the form ( |83| ) with A = 0, i.e. it is given by the tensor fields 
6|)-((79|). In order to write ^f lat in a convenient form, we introduce vector 



fields which are given in normal coordinates by 



v P = — P k d k d l r~ 1 



r Vp, 



Pp 



(r 2 P l -3x l P k x k ) e Q 2 



j J =e i *J j d k r- 1 = 



-3 i 

r Pj, 



1 



r~ 3 PA, 



P l j = -e ijk JjX k e Qi, 
P\ = -\ Ax ' e Si. 



(135) 
(136) 

(137) 



v* Q = -2 QV- 1 + -Q k d k d\ = r-%, 



-\r 2 Q l - \x l Q k x k e Q 2 , 
(138) 



where P l , J 1 , A, Q 1 are chosen such that the vector fields satisfy 



(C v P ) ab = * P °, (C Q vj) a0 = * 



yob 



ab 



. a b 



= * 



ab 



a b 



with f p b , *} b , 4^ b , ^ as given by We have on B a \ {i} 

(L Q v P ) a = 0, (L vj) a = 0, (L v A ) a =Q, (L v Q ) a = 0, 
and can thus write on B a /2 \ {i} 



^ A i 



(139) 



(140) 



= (C (v P + v,j + v A + vq)} 11 - - h ij h M (C (v P + vj + v A + v Q )) kl . 



' sing 



(141) 



Of the field <&° b g € W 1,P (S), p > 1, entering theorem 16 we assume that it 



can be written near i in the form 

$ ab = r s $" b + <j>" b 

reg ' Ireg 1 Zreg? 



(142) 
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where s < — 1 is some integer which will be fixed later on, <5?,! e „, $2reg are 

q(j,-s-i^ anc j XiXj&l — r 2 $ with 



smooth in B a and such that ^i^eg 
some $ G C°°(B a ). Then 

J reg — ^reg 1 ° * ^3 Y 2reg' 

with J 1 = r 2 Dj + s Xj G Qoo and J 1 = 0{ r - s ). 



(143) 



Using theorem 4.11. we obtain for the solutions of theorem |4.6| (where we 
can set by our present assumptions luq = 1, h = h') the following result. 



Corollary 4.12 With the tensor fields $° fc „ g; 9 iven by \l4A >, ^4A) re- 

spectively, let the vector field v a be such that 



*° 6 = *$na + *fea + ^v) a \ 



(144) 



satisfies D a * ab = in B a \ {i}. 



(i) If P a = in (14J-) and s = —3 in (141), ^ e vector field v a can be written 
in the form 

v l = r- 3 v\ + v\, with v\ € Qoo(B a ), v\ = 0(r 3 ), v\ € C°°{B a ). 



(ii) If J a =0, A = 0, Q a = in \i4iy and s = -5 in ([14%) , the vector field v a 
can be written in the form 

v* =r- 5 v\+vl with «ieeoo(B ), v\ = 0(r 5 ), v 2 G C°°(B a ). 

Proof: In both cases the vector field v a satisfies ~Lh,v a = — J sing ~ ^reg with 
J reg gi ven by (143) and J" in „ = Db®1i ng - By equation (14C) we have in case 
(i) J stng = {Lh(vj + v A + v Q )) a = (L h (vj + v A + v Q )) a , and in case (ii) 
Jting = (J->h,vp)) a = (LhVp)) a on B a / 2 \ {i} ■ The results now follow from 
equations (135)-( |l38| ), lemma [OO, and theorem 4.11. ■ 



We are in a position now to describe the behaviour of the scalar field ty a b 
near i. 



ab 



if P % P l +r- 2 h 2 



Lemma 4.13 The tensor field (\144\) satisfies 
in case (i) r 8 * ah * ah G E°°{B a ), 

in case (ii) r 8 ^ ah — ip+r ip R where ip R G C a (B a ) and ip 
with harmonic polynomial h 2 = ^r 2 (3 {P l ni) 2 — Pi P l ). 

Proof: For w 1 G Qoo(B a ) with Xi w l — r 2 w, w G C°°(B a ) we have 

(£ h (r s w)) ij =r s {{C h w) 13 - \sh ij w) +r s ^2^( 1 w J ). (145) 
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We set s = — 3 and w l = p l j +P\+Pq in case (i) and s = — 5 and w l — p l P in 
case (ii) and we write X\ v\ — r 2 v\ with v\ G C°°(B a ). Observing the equation 
above we get on B a / 2 \ {«} a representation 

i&i = r s II" + r s 2 K** + L ij 

with fields 

= (Co w) ij - H ft« fc w (Co w) kl ~^sw (s ij + 2 h ij (1 - i <5 fci ) 



+(£, iWl ) fci -^ s ^t) 1 + $i; eg , 

= 2 s ^ + v[ ] ), L ij = <S>g eg + (C h v 2 ) ij , 

which are in C co (B a / 2 ). Since a direct calculation gives Kij K l i = r 2 K with 
K G C°°(B a/2 ), we get 

* y >|/'-' = r 2s ~ 2 (2 //, , A '-' - K) + r 2s II,, H ij 

+r s - 2 2 I." + r s 2 //, ,/.'■' + /.,, L l \ 

from which we can immediately read off the desired result in case (i). In case (ii) 
it is obtained from our assumptions by a detailed calculation of r 2s_2 (2 Hij K v — 
K) + r 2s Hij H ij . U 

Combining the results above and observing the conformal invariancc of the 



equations involved, we obtain the following detailed version of theorem 1.2. We 



use here the notation of theorem 4.6 



Theorem 4.14 Assume that the metric h is smooth and \l/ ab is the solution of 
the momentum constraint determined in theorem \j.<\. If 



ft) ^ing = *J + *A + *Q ~ I h " b h cd (*J + *A + *Q ) m B a/2, 
(ii) $? eg = T-^teg + Keg ^th ^reg^teg € C°°{B a ) Such that 

0(r 2 ), and x a x a $^ es = ?" 2 * with some * <= C°°(B a ), 
then ^> ab satisfies condition luty. 



<•<) 



A On Holder functions 

In this section we want to proof an estimate concerning Holder continuous func- 
tions. 
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Let B be an open ball in R™, n > 1, centered at the origin. Suppose 
/ G C k (U) for some k > and m is a non-negative integer with m < k. Then 
we can write 

f= E i^/(0)a^ + m Al-tr 1 E j,d> 3 f(tx)xl 3 dt 

\P\<m P ' J ° |/3|=m " 



= E ^/(o)^+m | (l-o™- 1 E j { {^f{tx)-dPm)7fidt, 

\P\<m ' |/3|=m 

where the first line is a standard form of Taylor's formula and the second line a 
slight modification thereof. We denote by T m (f) the Taylor polynomial of order 
m and by R m {f) the modified remainder, i.e. the first and the second term of 
the second line respectively. 

Lemma A.l Suppose f e C m > a {U). Then f - T m (f) e C m < Q ([T) and we have 
for [3 € |/3| < m 7 

i5 /3 (/-T m (/))(x)i<|xr"+«-i' 3 i e V 7+/3 071 ^ (146) 

| 7 |=m-|/3| 7 ' 

where the constants c$ denote the Holder coefficients satisfying \d s f (x)— d s f(0)\ < 
c s \x\ a inU for S e Nft, |<5| = m. 

Proof: Applying the modified Taylor formula to / and then to its deriva- 
tives, we get 

d f3 {f - T m {f)) = T m _ m (dPf) + R m - m {d?f) - d^T m (f). 
We show that 

T m - W (d^f)-d^T m (f) = 0. (147) 

To prove this equation we use induction on n. For n = 1 the result follows by 
a direct calculation. To perform the induction step we assume n > 2 and show 
that the statement for n — 1 implies that for n. 

We write x = (x', x n ) for x € R" and (3 = (/?', f3 n ) for [3 e Nq etc. Then we 
find the equalities 

^T m (/) - ^>"(E T m _ 7B (9T»/) — r(x n )T») 



m-|/3|+/3„ 

= E r m _ l/3 ,|_ 7n (^'^/) 
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m-\i3\ 

Here the first line is a simple rewriting where we denote by T m _ 7n (<9 7 ™/) the 
Taylor polynomial of order m — j n of the function <9 7 ™ f(x' , 0) of n — 1 variables. 
In the second line the derivatives are taken and the induction hypothesis is 
used. The third line is obtained by redefining the index j n and using a similar 
rewriting as in the first line. 

With Jl47|) the estimate (|l46| ) follows immediately by estimating the integral 
defining R m -\f3\(d p f). ■ 



B An additional result 



In this section we prove a certain extension of theorem 3.11 



Theorem B.l Let u be a distribution satisfying Lu — f , where f £ E m ' a (B a ) 
and the coefficient of the elliptic operator L are in C m,a (B a ). Then 



r 3 ^2u k + u R £ E 

k=0 



m+2,a 



(B a ) 



(148) 



with u k e V k and u R £ C m+2 > a (B a ) 



Proof: We follow the proof of Theorem 3.11 using Schauder instead of LP 
estimates. Since f = fx + r/2 G E rn,a (B a ) we have 



f = r T m + fa with T m = ^ t k , 



k=0 

where T m is the Taylor polynomial of order m of ji and tk £ Vk- 
Consider the recurrence relation 

A(r 3 u )=rto, A(r 3 u k ) - r(t k - uj? \ 1 < k < m, 

which is obtained by defining k = 1, . . . , m, by 

k-l 

„3 \ " 



L(r*Y, u ]) = rC/(fe) ' 
3=0 



and defining and [/j m+1 ^ as in the proof of theorem 3.11. The equation 



for u and ( |148| ) then imply for ur equation (50) with s = 3. Since by our 
assumptions and lemma 3.E the right-hand side of this equation is in C m ' a (B a ), 
the interior Schauder estimates of theorem 2.6 imply that ur £ C m+2 (B a ). ■ 
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